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INTRODUCTION AND SUMMARY 


Radar backscattering from the sea is governed by two distinct physical 
processes. For viewing angles 0 near nadir the backscattering is primarily 
specxilar reflections frm sea waves with wavelengths much longer than the ra- 
diation wavelength X, The geometric-optics formulation of this large-scale 
scattering results in the normalized radar cross section c° being proportional 
to the probability density function P(n) of the large-scale surface normal n- 
For large nadir viewing angles the backscattering is principally Bragg scat- 
tering from sea waves with wavelengths comparable to X, and o° depends on the 
power spectrum W(ic,n) of these small-scale waves. The argiment ic denotes the 
vector sea wavenumber, and n denotes the normal of the underlying large-scale 
roughness. 

In this report we derive a two-scale scattering model that combines the 
two types of scattering in a maiuier consistent with energy conservation, that 
includes the effect of the tilting of the small-scale roughness by the large- 
scale roughness, and that accounts for the reduction of reflected power due to 
Bragg scattering. The special case of backscattering for which the transmitted 
polaiization equals the received polarization is then considered. An aniso- 
tropic large-scale surface of the type reported by Cox and Munk [l956] is used 
to specify P(n). In order to isolate the azimuthal variation of o® produced 
by the anisotronic P(n), we assume an isotropic small-scale spectrum. 

Computations of oo are compared with the AAFE RADSCAT data at 13.9 GHz for 
three wind speeds, U = 3, 6.5, and 15 m/s. Better agreement occurs for small 
and large 0 than for intermediace angles. The poor agreement in the midrange 
may be cause by the abrupt splitting of the sea spectrum into large- and small- 
scale components and by the somewhat arbitrary choice of 0.25 for the pertur- 
bation parameter. The 3 m/s data are considerably less than the computed values 
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at large nadir angles, and this nay be due to the friction velocity being less 
than the critical velocity required to significantly disturb the surface. At 
15 m/s the model displays an anisotropy that closely corresponds to the measured 
data. The anisotropy of the model decreases with decreasing wind speed* and 
the 3 m/s computations are essentially independent of the azimuth viewing angle 
In contrast, the anisotropy of the measurements does not change much with 
wind cpeed, and for light winds it is probably due to a directional small-scale 
spectr m, which is not considered in the present model. 

P?- rametric computations of a° at lU.6 GHz are also presented. The skewness 
of P(n) results in the ratio of the upwind to downwind o° oeing less than unity 
for small 0. This agrees with X-band measurements [Skolnik, 1970] except that 
the ratio for the measurements is smaller than that given by the model. The 
average <aO> of the upwind, crosswind, and downwind o® is also found. For small 
0, <a°> decreases with increasing wind speed, and at large angles the opposite 
is true. Tie average is least sensitive to wind speed variations near 15°. At 
30° the wind dependence is approximately for horizontal polarization and 
for vertical polarization when U i 20 m/s. For large angles the hori- 
zontally polarized <a°> has a stronger wind dependence than does the verticsdly 
polarized <a°>. Finally, the dependence of <o°> on 0 and on polarization di- 
minishes for strong winds. 



TWO-SCALE SCATTERING MODEL 


The mean sea surface defines the z = 0 plane of a x,y,z coordinate system, 
and the x axis points in the upwind diroction. All vectors are referenced to 
this system. Unit vectors are denoted by a caret, and vectors that in general 
do not have unit magnitude are denoted by an arrow. The axis vectors are de • 
noted by X, y, and z. The components of a vector are indicated by deleting 
the caret or arrow and inserting the superscript x, y, or z. The radar bore- 
sight direction is represented by pointing towards the sea surface I , and 

the viewing angle 0 is defined as the angle made by and z. 

The definition of scattering coefficient, which will be used throughout 
this report, is the ratio of the power density p(kg) to the time-averaged 
power incident onto the surface in question, where p(kg)dk2dkg is the time- 
averaged scattered power having a propagation vector iri the neighborhood dkgdkg 
of kg. This is *not the usual definition but is more compatible witn the k^ ,ks 
vector notation used herein. Note that the integrcd of the scattering coef- 
ficient over all ks equals the ratio of the total scattered power to the in- 
cident power. The scattering coefficient defined by Peake [l959] is in terms 
of power scattered per solid angle and is f.->Tind by multiplying our coefficient 
by Uirkl (a differential solid angle dfi = dk^dk^/k|). The normalized raiar 
cross section o^(kj^) is given by the product of Peake’s scattering coefficient 
and cos 9. Hence in terms of our scattering coefficient r(kj^,kg) for the sea 
surface, the normalized radar cross section is 

a^(k^) = cos^O r(ki,-kj) (l) 

A two-scale scattering model is used to compute r(kj[,kg). The sea surface 
is modeled by a small-scale surface Eg superimposed onto a large-scale surface 

The two-scale surface E is then the sum of E3 and E^. 


The rms height 



u 


variation on Zg is assumed small compared to the radiation wavelength A, and 
the rms slope variation is assumed small compared to unity. These two re- 
quirements are necessary for the application of perturbation theory [Rice, 
1951] in treating the radiaticn scattered by Zg. The radius of curvature 
at all points on Z^ is assumed much greater than A. This allows for dividing 
Z^ into finite surface elements AZ^ that have dimensions large relative to A 
and that are nearly flat in the respect that the variation of the normal to 
AZ^ is small compared to the mean normal n of AZ^. The electric field on the 
two-scaile surface element AZ associated with AZ£ is approximated by the field 
that would be present on the infinte plane normal to n* Zg being superimposed 
on the plane. Furthermore, the fields on adjacent elements are assumed Tin- 
correlated. This results in the scattering coefficient of AZ being equal to 
the scattering coefficient r(kj^,k 3 ,n) of the tilted small-scale surface Zg, 
and the total scattered power is the sum of the power scattered by the indi- 
vidual elements. Multiple scattering is not considered, and all of the scat- 
tered power is assumed to escape from the surface. 

The power density p(koJ for the two-s tie surface is found by summing 
over all surface elements that are illuminated by the incident radiation. 

p(ks) = I]r(ki,ks,nni) Api(ki,rini) (2) 

m 

where is the mean normal for the m*^^ element and Ap^(k^,fijj^) is the power 
incident onto the element. Divridlng (2) by the total incident power gives 
the scattering coefficient for the two-scale surface. 

r(ki,kg) = (A cos e)”^ 2Zr'(ki,kg,%) (3) 

m 

where A is the area of the mean two-scale surface that is subtended by the in- 
cident plane wave and AAj^^ is the area of AZ^ for the element. The 



assumption that the dimension of AE^ is small relative to means that the 
variation in n from one element to the next is small, and hence the above 
summation can be replaced by the integral 
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r(ki,kg) = (A cos 0) ^ /dn* /dn^ (-ki«n) 5(k£,fi) r(k^,kg,n) (U) 

where C(ki >n)dn*dn^ is the area of that has a normal in the neighborhood 
dr\ dn of n and that is illuminated by the incident radiation. The region 
of integration is (n^)^ + (n^)^ ~ !• 

Let H(n)dn^dn^ denote the area of that has a norms'*, in the neighbor- 
hood dn^dn^ of n and that is subtended by the inc.* lent wave. The probability 
density function P(n) of the large-scale siirface normal is then defined by 

P(n) = n- E(n)/A (5) 

Let l(kj,n) be the fraction of E(n)dn^dn^ that is not shadowed by a remote 
portion of the surface from the incident radiation. The illuminated area den- 
sity function 5(k^,n) is then given by the product of H(n) and l{k^,n). 

5(k^,n) = l(kj^,n) P(n) A/n^ (6) 

The illumination function l(kj^,n) is readily found by assuming that it is in- 
dependent of n except through the unit step function uC-kj^^n), which accounts 
for the situation in which the angle between n and exceeds tt/ 2, totally 
ruling out the possibility of illumination. 

I(ki,n) - x(ki) u(-ki*n) (7) 

The function x(^i) is detei^^^ned as follows. The ratio of the power incident 
onto area C(kj^,fi)dr)^dny to the total incident power is given by 


r(ki,n) = (a cos 0)“^ (-ki*n) C(ki ,n)dn^dn^ 


( 8 ) 
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The integral of r(k^,n) over all n is unity, and solving for one obtains 


x(&i) = cos e//dn* /dn^ G(kj^,n) 

(9) 

G(ki,n) = (-k^*ri) u(-kj^*n) P(n)/n^ 

(10) 


For situations in which shadewing is insignificant, such as smooth seas or 
small 0, x(^i) approaches unity* Substituting (6) and (7) into (4) yields 

r(k^,kg) = xi^i) sec 0 /dn^ /dn^ G(k 5 ^,n) r(k^,ks,n) (11) 

In the Appendix r(kj^,k 3 ,n) is derived as the sum of two scattering coef- 

ficients, one associated with the incoherent scattered power and the other with 
the coherent reflected power. The coherent coefficient given by (A25) has the 
form of a Dirac delta function, and when it is substituted into (ll) yields 

r^(kj^,kg) = x(ki) sec 9 P(no) Y(k^ ,fio)/Ukg (12) 

% = - kil (13) 

where the power reflection coefficient y(k^ is given by (A2U). 

We now consider the special case of backscattcring, i.e., kg = -k^, for 
which the transmitted polarization equals the received polarization. Also the 
roughness on the small-scale surface Zg is assumed isotropic. Under these 
conditions (12) becomes 

^ x(ki) sec^e P(-k.) y(ki,-ki) (U+) 

where Y(ki,-ki) is given by (A32). The incoherent two-scale backscattering 
coefficient is given by 


r^(ki,-ki) = x(ki) sec 0 /dn^ /dn^ G(kj_,n) F^(kj ,-kj^ ,.l) 


(15) 



where r^(kj^,-kj^,n) is given by (A3l). Summing the coherent and incoherent co- 
efficients gives the two-scale backs catte ring coefficient. 


r(ki,-fci) = r^(ki,-ki) + r^(%,-k^) 


(16) 
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SPECIFICATION OF SEA SURFACE ROUGHNESS 


The two-scale scattering model requires as inputs two distributions that 
characterize the sea surface roughness. These are the small-scale roughness 
power spectrum W(k), ic being the vector wavenumber, and the probability den- 
sity function F(n) of the large-scale surface normal n- For the purpose of 
isolating the effect of an anisotropic large-scale surface, we assume an iso- 
tropic small-scede spectrum W(k) that depends only on k = |ic| . Pierson and 
Ctacy^s [1973] empirical sea spectrum S(ic) (Eqs. 2.5 - 2.9 in their report) 
is used to specify W(k). The amplitude of S(k) is a function of the friction 
velocity u^, which in turn is a fianction of the wind speed, anemometer height, 
and air-sea temperature difference [Cardone, 1969]* 

The spectrum S(k) is divided into a large-scale spectrum and a small- 

scale spectrum Sg(K). 


s,(0 = 


S(k) ic < k. 


>C > K, 


,(<) = 


0 

s(<) 




K > K, 


and W(k') is related to S_(k) by 


W(ic) = i.2h) S_(*c)/k 


(17) 


(13) 


(19) 


The value of the cutoff wavenumber is found by assigning a value to the 
small-scale perturbation parameter kc, where k is the radiation wavenumber and 
r is the nns height variation on the small-scale surface. Integrating over 
the small-scale spectrum gives 

00 

- (^/2) JdK K W(k) ( 20 ) 

0 



or in teras of S(k) 


= fdK S(k ) (?1 ) 

<c 

iTie above int-'->r,rul can be evaluated in closed form, ^ being expressed as a 
function of and u^^. The inverse of this function gives in terms of i; 
and u,j. Perturbation theory requires that kc be small in comparison to unity 
and setting equal to zero results in the two-scale model degenerating to 
geometric optics. We use an intermediate value of 0.25 for the radar cross 
section computations, and the values of ic^ appear in the next section. 

The specification of the large-scale slope density P(n) is based on Cox 
and Munk's [l956j measurements of the sun glitter on rough seas. Their dat.n 
were reduced in terms of the probability density function P(Zy,Z£.) of the 
up /downwind and crosswind surface slopes, Z^ and Z^. The relationship be- 
tween P(n) amd P(Z^,Zg) is 


P(n) = P(Zu,Zc) 

(22) 

Zy = -n*/n^ 

(23) 

Z, = -n^/n=^ 

(2U^ 


where is the Jacobian relating the Zy,Z(, coordinates to the co- 

ordinates. 

The sun glitter data were fitted to a two-dimensional Gram-Charlier serie 
and P{Zy,Z^) was found to be close to Gaussian with some up/downwind skewness 
that increased with wind speed. The most probabie slope for high winds was 
about Z,j = -tan i°. 'The data also showed a peakedness, barely above the exper 
Imental error, such that the probability of very large and very small slopes 
was greater than Gaussian. These properties are explicated in the following 
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expressions jiiven by Cox and Munk: 

P(Zy,Z^) = [l + T(u,v)] exp[-ig(p2 + v2)] /(2ir<zJ>^ <zj>^) (25) 

T(u,v) * Cjp(v^-l) + CjCu^-Sb) + C3 {v**-6v^+3) 

c^(v2-l)(|i^-l) + C5 (u**-6u^+3) (26) 

w = (-7) 

V = '?8) 


2 2 

where <Z^> and <Z^> are the up/downwind and crosswind slope variances. The skew- 
ness coefficients are functions of the wind speed U (in m/s) 


Cj = -(0,01 - 0.0086 u )/2 (29) 

C2 = -(O.OU - 0.033U)/6 (30) 

and the peakedness coefficients aa*e constants. 

C3 = O.UO/2U (31) 

c^ = 0.12/U (32) 

C5 = O.23/2U (33) 


The sun p;litter from large and infrequent slopes was masked by a back- 
ground of sunlight scattered by submerged particles and reflected skylight, 
and the vaiues for the slope variances reported by Cox and Kunk represent a 
lower bound [Wentz, 1975], To correct for this .nd to filter out the slope 
contribution of the smali-scede roughness, we multiply the Cox-Munk up/downwind 
and crosswind variances by the rat’o of the total variance calculated from the 
Pierson-Stacy large-scale spectrum to the Cox-Munk total variance, ''he var- 
iances that are used in (25) sire then 
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<z?> 

= (<z^> /<z^> 

) 

<7?> 

U 

p-s c« 

-m 

u c-m 

<z|> 

= (<z^> /<z^> 

) 

<'2> 


p-s C- 

-m 

^ C-m 


where the values of the variances reported hy Cox and Mimk are 

<4>c-m = 3 .i6x10"3u 

<Zc> = 0.003 + 1.92x10* 
c-m 

and the total variance is 

2 2 2 
c~m n c-m c c-m 

Tlae Pierson-Stacy total variance is calculated from 


(3h) 

(35) 

(36) 

(37) 

(38) 


<Z^> g = fdK S(ic) 
o 


(39) 


Although In the strictest sense (39) only applies to a Gaussian surface, the 
deviation from Gaussir . indicated by the skewness and peakedness is slight and 
probably does not cause s- mi'^icant error in the variance calculation. Note 
that the Pierson-Stacy variances for the entire ser spectrum (found by inte- 
grating from 0 to ») are about twice as large as the Cox-Munk variances. The 
values of the total variances used in the radar cross section computations 

(i.e. the sum of (3^) and (35), which equals <Z^> ) appear in the next sec- 

p-s 

ticn. The error in skewness and peakedness due to the lack of sun glitter data 


for large slopes is not considered. 
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COHPUTATIOHS OF THE NORMALIZED RADAR CROSS SECTION 

llie backscatterizig model described in the preTiotis two sections requires 
the following inputs: (l) radiation frequency f, (2) permittivity e of sea- 

water. (3) wind speed U. (U) nadir viewing angle 9. and (^) azimuth vieinng 
angle 4* frequencies are considered: 13 >9 GHz for the conparisons with 

the AAFE RADSCAT data and li|.6 (Sz for the parametric computations. The per- 
mittivities for these two frequencies are shown in Tables 1 and 2 and are cal- 
culated from expressions given by Porter and Wentz [l97l] for a seawater tem- 
perature and salinity of 28A<>K and 33 °/oo. 

The sea surface rou^mess distributions W(k) and P(n) discussed in the pre- 
vious section depend on the wind speed U and friction velocity u,. Cardone's 
[ 1969 ] eiqpressions for a neutrally stratified atmosphere and for an anemometer 
hei^t of 19*9 inters are used to calculate u,^ as a function of U. Values 
for u« along with the cutoff wavenumber and the large-scale slope variance 
<Z^>, which are referred to in the preceding section, also appear in Tables 1 
and 2. 


Table 1. Inputs for the AAFE RADSCAT Cosparisons 
Frequency = 13-9 GHz 
Permittivity = UO.l - 39- 3i 

Wind Speed Friction Velocity Cutoff Wavenumber Slope Variance 
U u* Kc <Z^> 

(m/s) (cm/s) (cm~l) 


3.0 

6.5 

15.0 


10. U 

21.5 

61.9 


0.79 

0.97 

1.63 


0.018 

0.032 

0.097 



13 

Table 2. Inputs for the Paranetric Computations 
Frequency = lU.6 GHz 
Permittivity = 38. U - 39. Oi 


Wind Speed 

Friction Velocity 

Cutoff Vavenuniber 

Slope Variance 

U 

«a 

•'c 

<Z2> 

(m/s) 

(cm/s) 

(cm“l ) 


5 

l6.0 

0.93 

0.026 

10 

37.2 

1.27 

0.054 

15 

61.9 

1.71 

0.099 

20 

91.1 

2.28 

0.175 

25 

122.4 

2.93 

0.293 

30 

158.5 

3.75 

0.492 


The azimuth viewing angle ^ is the angle made by the x axis, which points 
upwind, and the projection of the bores i^t vector kj onto the z = 0 plane, 
which is the mean sea surface. The boresight direction is then specified by 

ii^ - (cos ^ sin 6, sin f sin 6, -cos 6) (UO) 

and f = 0^, 90°, and l80° refer to measurements looking upwind, crosswind, and 
downwind, respectively. The backscattering model has no crosswind asymmetry, 
and hence the ccmiputations for ^ and are identical. The incident polari- 
zation vector Pj, which appears in the Appendix, is calculated from 

k.x£/|k{xz| horizontal polarization 

ki»«(£.^xz)/|fcix£| vertical polarization 

Computations are done for the following directions: 

0 = 0°, 10°, 15*^, 20°, 25°, 30°, 40°, 50°, 60 ° 

♦ = 0 ®, 90 °, 180 ° 


In addition, ccmiputations for 13.9 GHz and for 0 » 30° are done for ^ ranging 
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from 0® to 180® in 10® steps. In the figures these cooputation points are con- 
nected by strai^t lines. 

The model is compared vith the AAFE RAI^SCAT data at 13.9 GHz in Figures 1 
tbrou^ 5. The horizontally polarized normalized radar cross section is 
plotted versus the nadir viewing angle 6 in Figures 1 through 3 for upwind, 
crossvind, and downwind observations. The short-dashed, long-dashed, and solid 
curves represent the computations for 3, 6.5, and 15 m/s winds, respectively. 
The measurements for these three wind speeds are indicated by squares, circles, 
and stars, respectively. The agreement between the model and the measurements 
is better for small and large values of 6 than for intermediate values. At 
small nadir angles the backscattering is primarily specular reflections from 
the large-scale surface, whereas at large angles small-scale Bragg scattering 
dominates. The poor agreement in the transitional region may be due to the 
abrupt splitting of the sea spectrum into large- and small -scsde components and 
to the somewhat arbitrary choice of 0.25 for the pertxirbation parameter. 

The model displays a stronger wind dependence near nadir than do the meas- 
>jrements. At large angles, for which the opposite is true, the agreement at 
15 m/s is fair, but the 3-5 m/s mesisurements are considerably less than the 
ccmputations . Experiments in wind-water tunnels show a sudden increase in 
wave height at a critical friction velocity near 12 cm/s [iterson suid Stacy, 
1973^ , and the expressions that are used to specify the sea spectrum are valid 
only for u, greater than this critical velocity. The low values of the 3.5 m/s 
data at large angles possibly indicate that the friction velocity during these 
mesisurements was less than the critical velocity. This possibility is sup- 
ported by the u,| calculations, which give a VEdue of 10. U cm/s for a 3.5 m/s 
wind speed. 

In Figures U and 5, o® at 6 = 30® is plotted versus the azimuth viewing 
angle ^ for horizontal and verticsd polarizations. The same convention as 
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Fig. 3. Dovmrlnd normalized radco* cross section for horizontal polarization 
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Pig. 1». Azinutb variation of the normalized radar cross section for horizontal polarization 
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Fig. 3. Azimuth vmria'' ion of the normalized radar cross section for vertical polarization 
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described above is used to identify the computations and measurements at the 
three wind speeds. A blanket 6 db has been subtracted from all computations 
in order to align them with the measurements • The spacing between the three 
wind speed curves agrees ve*i with the observations. The 15 m/s curve shows 
a significant anisotropy that closely corresponds to the measured data. This 
dependence of the model on ^ is due to reflections from the anisotropic large- 
scale surface. The ^ dependence decreases with decreasing wind speed, and the 
3 m/s curve is essentially flat. In contrast, the anisotropy shown by the 
measurements does not change much with wind speed, and at low wind speeds it 
is probably due to a directional small-scale spectrum. 

Parametric computations for lU.6 GHz are presented in Figures 6 through 
11. The upwind-crosswind and upwind-downwind ratios of a® are plotted versus 
9 in Figures 6 and 7 for horizontal and vertical polarizations. Computations 
for three wind speeds, 5» 10, and 20 m/s, are shown by the short-dashed, long- 
dashed, and solid curves, respectively. X-band measurements, which are re- 
ported in the Radar Handbook [Skolnik, 1970] , are indicated by crosses euid rep- 
resent an average over six days during which the median wind speed was 5 m/s. 
The measured upwind-downwind ratio drops sharply at 8 = 10^. In comparison, 
the model shows a dip but not as extreme because of the small skewness in the 
large-scale slope probability. For the larger angles the upwind-downwind ratio 
becomes greater than unity for both the computations and observations, although 
the 5 itt/s curve is essentially at 0 db. The discrepancy at large angles is 
probably due in part to neglecting the small-scale anisotropy, which seems to 
be an important factor for light winds. 

The average <o®> of over is plotted versus 0 in Figures 8 and 9 for 
horizontal and vertical polarizations, 

<aO> . ^5(o°p + 


( 42 ) 
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Pig. 6. Upwlnd-crosBvind and upwind-downwind ratios for horiaontal polarization 
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Fig. 7. Upwind-crosswlnd and upwind-downwind o® ratios for vertical polarization 
















Six wind speeds, ranging from 5 to 30 m/s in 5 ra/s steps, are shown. The 
curves cross each other between 10® and 30o, bu+ ist considering the 5: 10, 
and 15 m/s curves, the crossover region is more narrow, being between 10'-' ana 
15^* The curves flatten out considerably with increasing wind speed, and 
at 30 m/s drops only about 8 db from 0® to 60 ^ and is nearly independent of 
polarization. The polarization independence is due to the dominance of large 
scaJ.e reflections for strong winds. 

The same computations that appear in Figures 6 and 9 appear again in 
Figures 10 aiiu 11 except that they are plotted versus the log of wind speed 
rather than 6. Curves for ten nadir viewing angles are shown. Near nadir 
<o^> decreases with increasing wind speedy and for the larger angles the op- 
posite is true. The curve for 0 = 15° least sensitive to wind speed vari- 
ations. At the ^arger angles the horizontally polarized <a®> has a stronger 
dependence on wind speed than does the vertically polarized <a®>. At 30® 
the wind dependence is about the same, being for horizontal polarization 
and for vertical polarization when U £ 20 m/s. 
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Pig. 11. The vertically polarized <<r®> versus the log of the wind speed 
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APPENDIX: SMALL-SCALE SCATTERING COEFFICIENT 


According to Rice's [l95l] perturbation theory, a plane wave incident 
onto a random, slightly rough surface E produces an incoherent scattered 
field and a coherent reflected field. The scattered field is represented by 
a set of plane waves, and the electric field E(kg) of one such plane wave 
having a propagation vector kg is (suppressing the time dependence) 

E(ks) = [{Pi-Hi)(BhhHs + Vs) + (Pi*Vi)(BvhHs B^Vs)]exp[ik(£s*?)] (Al) 

where k is the radiation wavenumber, r is the position vector, end Pj is the 
polarization vector of the incident field. The inc dent horizontal and ver- 
tical polarization vectors. Hi and Vi, that €tre referenced to the normal n to 
the mean surface are given by 

Hi = ki*n/|kix»^j (A2) 

* k^xH- (A3) 

where k^ is the propagation vector of the incident plane wave. Th-i scattered 
polvrization vectors, Hg and Vg, are given by (A2) and (A3) with the subscript 
s replacing the subscript i . Peake and Barrick [ 1967 ] derived the scattering 
terms m = h or v and n = h or v, to first order in the pertvurbation pa- 

rameter k?, where c equals the rms surface hei^t variance. 

Bmn - -2k(-ki*n)onnI*('^) (aU) 

where N(ic) is the coefficient of the roughness spectral component having the 
vector wavenumber ic. 


< = kj(ks-ki) - [(ks-ki)*n]n I 


(A5) 
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The bistatic matrix elements Oqq given by Peake and Barrick in terms of 
the surface permittivity e and the angles 6x, 6g, and *g. Thes' angles are 
related to kj^, k^, and n by the following equations: 


= arccos(-ki*n) 

(A6) 

= arccos(kg*n) 

(A7) 

= arccos(Hi*Hg) 

(A8) 


We use the convention that the first a subscript refers to the incident polari- 
zation and the second subscript refers to the scattered polarization. Peake and 
Barrick used the opposite convention. Also our Hg is the negative of that de- 
fined by them. Accordingly the following are the eq>propriate substitutions: 


®hh ' 

(A9) 

II 

(AlO) 

a . = -a ' 
vh hv 

(All) 

a = a* 

w w 

(A12) 


where the are those appeairing in Peake and Barrick [l9o7] . 

A ^ 

Pg denote the polarization vector of the receiver. The Pg polari- 
zation component of the time-averaged power of plane wave kg divided by the 
time-averaged incident power is given by (* denotes complex conjugate) 

Y(ki,ks.n) = (|E(ks)-Ps|^)(kg*fi)/(-ki-n) (A13) 


where the angle brackets denote average over time. The surface roughness is 
assumed to e^erience random flucuations in time such that the time-averaged 
electric fields of the scattered plane waves are uncorrelated, and as a result 
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the powers of tne individual scattered plane waves are additive. The assun^)- 
tion of random roughness also implies 

<|N(i?)|2> = kJ W(iT) (AlU) 

where W(<) is the roughness spectrum defined by Rice and ic^ is the wavenumber 
of the fundamental roughness spectral coiDponent. 

The assumption that k is much larg'^r is made, and as a result 

the scattered waves are close together in kg -space, with the ::patjing being 

~ (AI 5 ) 

The kg-space distribution of scattered power is approximated by a continuous 
distribution for ^ich the power of the discrete plane waves is evenly spread 
over the spacings given by (A15). The scattering coefficient for the inco- 
herent power is then given by 

r^(kj,kg,n) = (AI 6 ) 


Combining the above equations yields 


rj^(kj,kg,n) - k**(-k^ 'nHitg’n)^ w(ic) |T|2/kg 


(AI7) 




The coherent reflected field is represented by a single plane wave prop- 
agating in the specular direction kj.. 

kj. = ki + 2 (-k.*n)n (AI9) 

The reflected electric field is given by 


E(kp) = [(Pi*Hj^)Rh ”r (Pi‘Vi)Ry Vr]exp[ik(k^*?)] 


(A20) 
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where the reflected horizontal and vertical polarization vectors, and V^., 
are given by (A2) and (A3) with the subscript r replacing the subscript i. 

Wu and Fung [l972] expressed the horizontal and vertical polarization i*eflec- 
tion coefficients, and Ry, in the form 

Rm = Pm^0i) - <3m) (A21) 

oo oo 

Qjj, = (k cos 0^/2) /du /dw W(l?) Fj„(u,w) (A22) 

_00 ^00 

where m = h or v, are the Fresnel reflection coefficients, and the 

functions Fjjj(u,w) are given by Wu and Fung and have an implicit dependence 
on the surface permittivity e, k, and 0^. The roughness vector wavenumber ie 
in this case is given by the two-dimensional vector 

ic = (u - k sin 0^, w) (A23) 

Note that 0^^^ is of order (k^)^. 

The Pg polarization component o'f the time-averaged reflected power di- 
vided by the time-averaged incident power is 

y(ki,n) = jE(kr)*Psl^ 

A ^ 

The scattering coefficient for the coherent power equals zero for all kg ^ k^, 

and its integral over all kg equals Y(l^i»b). It thus has the form 

ro(ki,kg,n) = 6(k* - kj) 6(k^ - kj) (A25) 

where 6 denotes the Dirac delta function. The total scattering coefficient 
for the random, slightly rough surface Zg is then 


r(kj,kg,n) = rx(kj^,kg,n) + rQ(ki,kg,n) 


(A26) 
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A 

We now consider the special case of backscattering, i.e., kg = -kj^. 

The bistatic matrix elements for backscattering are 

Mjjjj = (l - e)/fcos 0j + (e - sin (A27) 

«hv = 0 (A28) 

«vh = ° ^^9) 

cxyy = (e - l)[(e - l)sin^9j^ + e]/[e cos 0^^ + (e - sin^0j^ (A30) 

We also require that the transmitted polarization equals the received polar- 
ization, i.e., Pj = Pg, and that the roughness spectrum W(ic) be isotropic and 
depends only on . Under these conditions the incoherent backscattering 
coefficient becomes 

r^(ki,-ki,n) » k‘*cos30i W(2k sin ©i)! - |Pj-Vi|2a^|2/(_kZ) (a31) 

The scattering coefficients are integrated over n to obtain the two-scale 
scattering coefficients. The two-scale coherent backscattering coefficient 
contsdns the term y(kj[,-kj^). Under the above stated conditions (A2U) takes 
the following form when -kj^ is substituted for q: 


Y(ki,-ki) = |p(0)|2 (1 - 2 Re Q) (A32) 

00 

Q » (wk/2) /du u W(u) [2ke^ ♦ 2b - 2c + u2(c - b)(u2 + bc)“^] (A33) 

0 * 

b = (k2 - u2)^ (A3*t) 

c = (ek^ - u^)^ (A35) 


where b is either positive real or negative imaginary. The Fresnel reflection 
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coefficient for normal incidence is 

p(0) * (1 - e^)/(l + e^) (Aj6) 

where in this case the sign is eurbitrary. Note that in (A32) the |qP term, 
which is of order (kc)**, has been dropped. 

Wentz [l97*»] proves that the above formulation satisfies energy conserva- 
tion to second order in k(. The proof entails taking the sum of the scattering 
coefficients given by (A26) for two orthogonal scattered polarizations. The 

A 

integral of this sxim over all kg is then shown to equal unity for a perfect 
conductor. 
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